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Abstract. We explore the intimate relationship between quantum lithography, 
Heisenberg-limited parameter estimation and the rate of dynamical evolution of 
quantum states. We show how both the enhanced accuracy in measurements 
and the increased resolution in quantum lithography follow from the use of 
entanglement. Mathematically, the hyperresolution of quantum lithography 
appears naturally in the derivation of Heisenberg-limited parameter estimation. 
We also review recent experiments offering a proof of principle of quantum 
lithography, and we address the question of state preparation and the fabrication 
of suitable photoresists. 
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An important branch of quantum mechanics is parameter estimation. Heisenberg's 
uncertainty principle seems to prevent us from determining a physical parameter such 
as a phase with infinite precision, and it is therefore important to understand what 
are the limits of the estimation process. It is also a branch of physics that particularly 
interested H.A. Haus (e.g., Haus, 1995). 

In this paper we consider the link between phase estimation and the rate of 
dynamical evolution of quantum states. This will lead us to the concept of quantum 
lithography and the question of the role of entanglement. Finally, we will review some 
experiments that show the viability of quantum lithography (at least in principle), 
and briefly consider the generation of the necessary optical quantum states. 

1. Parameter estimation 

Consider a physical process that induces a relative phase ip in the state of a system. For 
example, one can send light through a medium with an unknown index of refraction 
rir to induce a phase shift (p = 2Tr L{nr — 1)/A, with A the wavelength of the light 
and L the length of the medium. This phase shift can be measured relative to a 
reference beam in a Mach-Zehnder interferometer. When we assume that L and A are 
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known with very high precision, we can infer the index of refraction with an accuracy 
that is proportional to the error in the phase. Many precision measurements can be 
reformulated in terms of a phase measurement, with the search for gravitational waves 
as one of the most urgent cases (Fritschel, 1998). It is therefore important to have 
definite bounds on the error in the phase, given a particular set of resources. We 
should realize, however, that there are two distinct questions we can ask about the 
uncertainty in the phase. 

First, there is the relative phase difference dip that renders two quantum states 
IV'(O)) and \'ip{S(p)) distinguishable, i.e., 

{mm^)) = o- (1) 

Here, we assume that the dependency on S(p actually does render \ip(6ip)) orthogonal to 
|'0(O)) . We can then minimize Sip to satisfy Eq. ^ . This depends on the characteristics 
of IV'(O)) and liplSip)). Thus, Sip is a measure of the dynamical rate of the (unitary) 
evolution. 

A second question we can ask, is how well one can measure ip given a state \ip{ip)). 
That is, we want an expression for the error Aip in the phase. Since there is no phase 
operator, we usually measure a suitably chosen (Hermitian) operator X such that 

= . (2) 

\d{X)/d^\ 

Mathematically, these two problems are very similar, since measuring \ip{ip)) typically 
involves interference with some other state \ip{ip')). Choosing situations where 
Sip = Aip allows us to derive bounds on the fundamental quantum limits of parameter 
estimation (Aip) by calculating the maximum rate of dynamical evolution {Sip). 

It was recognized early on by Mandelstam and Tamm in 1945 that the uncertainty 
relation between time and energy for a given frequency uj can be used to prove the 
following inequality (Mandelstam and Tamm, 1945): 

Aw = ujAt > — - — , (3) 
^ - 2 AE ' ^ ' 

where t is the time it takes to evolve from the initial state to the next orthogonal 

state, hu! is the energy difference between the two orthogonal states and we have 

used Sip = Aip. Furthermore, AE is the uncertainty in energy. This inequality can 

be interpreted as follows: the minimum phase difference Aip that can be detected 

is inversely proportional to the normalized energy spread of the state l^p) (with 

{AEf = (Vl-B^ji/)) - {iIj\E\^)^). The normalization is given by hw. The factor 7r/2 

has a geometric interpretation that is highly relevant for this paper: The minimum 

distance between two points that can be distinguished by a wave is given by the 

distance between a crest and the adjacent trough of the wave. This (normalized) 

distance of tt/2 is also known as the Rayleigh limit. 

However, it became clear that the Mandelstam- Tamm (MT) inequality is not 

applicable to situations involving certain pathological states (Shapiro et al, 1989; 

Braunstein et al, 1992; Uffink, 1993). When coherent (classical) states with AE = 

hio^J {n) are used, the MT-inequality yields the Poissonian error 

vr 1 
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where (n) is the average number of quanta in the coherent state. The indication was 
that phase estimation could not improve beyond the so-caUed Heisenberg hmit (Ou, 
1996): 

A^>^. (5) 

But there exist states with finite average energy E and unbounded AE. According 
to the MT-inequahty, such states could have a precision that is proportional to, for 
example, 

If there are states for which the MT-inequality does not give a saturated bound, 
then the next question is whether an additional bound can be derived. In 1998, 
Margolus and Levitin did just that (Margolus and Levitin, 1998). They showed that 

At>--, (6) 
which for narrow-band schemes around frequency lu with E = hiij{n) yields 

A^>|y\. (7) 

This inequality therefore restricts the precision in a phase measurement by the average 
energy^ rather than the energy spread. The pathological states mentioned above turn 
out to have sufficiently small average energies so that Eq. ISJ is still satisfied. Recent 
examples of applied phase estimation include frequency measurements (Bollinger, 
1996), the quantum gyroscope (Dowling, 1998), quantum positioning and clock 
synchronization (Giovannetti et al., 2001), and length and weak force measurements 
using coherent optical states (Ralph, 2002; Munro et al., 2002). 

2. The role of entanglement in parameter estimation 

Suppose we want to estimate a parameter ip by measuring a particular observable. 
Construct the state 

|(p) = i=(|0)+e^^|l)) , (8) 

where the basis { |0) , 1 1) } is chosen in some convenient way determined by the physics of 
estimating (/? with the choice of observable as fJa; = |0)(1|-|-|1)(0|. Then the expectation 
value of Ux is given by {if\ax\'fi) = cosif. When we repeat this experiment N times, 
we obtain 

N 

(af ) =1 (^1 . . . e . . . W)m = iVcos^ . (9) 

fe=i 

Furthermore, we know that ct^ = 1, and the variance of Gx given N samples is readily 
computed to be {AdxY — ^(1 ~ i^) = A^sin^ Lp. According to estimation theory 

\d{'^x)lM 

The standard variance in the parameter p after A^ trials is thus given by 

A^, = ^""^ = 1 . (11) 

In other words, the uncertainty in the phase is inversely proportional to the square 
root of the number of trials. This is the classical Poissonian error in the phase. 
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With the help of quantum entanglement we can achieve the Heisenberg limit of 
1/N. Consider an entangled input state on N systems: 

I^at) = -^(|0,...,0)i...^ + e^^'^|l,...,l)i...^) . (12) 

The relative phase e*^''' can be obtained by a unitary evolution |1) e*'^|l) and |0) 
|0), thus yielding the required factor. When we suggestively write |0) = |0, . . . , 0)i...Ar 
and |1) = |1, . . . , 1)1. ..AT, then the state becomes 

|^^) = -i=(|0) + e^^^|l)) . (13) 

This is mathematically equivalent to a single (nonlocal!) system with a relative phase 
shift of Nip. In order to measure this phase, we need to measure the (nonlocal) 
observable Sat: 

EAr = |0>(l| + |l)(0| . (14) 

This yields 

{iPn\^n\(Pn) = cosiV.^ . (15) 
As before, we obtain 

ASjv ^ 1 sinNifi ^ 1 
|d(SA,)/d(^| NsinNip N ' ^ ' 



Here we see that the precision in ip is increased by a factor over the standard 
noise limit when we exploit quantum entanglement. 

When the loss of subsystems is considered (for example, when n out of N photons 
fail to arrive), the entangled state in Eq. H12(l becomes separable and mixed: 

PN-n = i|0)i...JV-„(0| + i|l)i...Ar_„(l| • (17) 

Not only is it separable, there is also no information about ip in this state. This state 
can therefore not be used for the parameter estimation at all. 

One way to circumvent this practical difficulty is to use separable states and non- 
local observables. When N — n out of N systems arrive, the experimenter chooses 
to measure Ejv_„ instead of S^r. This way, the parameter (p is estimated with an 
enhanced precision proportional to 1/{N — n). The price to pay is that often the 
measurement of Ejv-n will fail because it does not span the complete state space. 
It must also be noted that this trick depends on the physical implementation, and 
cannot always be applied. 

3. Quantum lithography 

One particular physical implementation of this process involves light, and this will lead 
us to the main topic of this paper. Suppose that the state in Eq. (|13|l is physically 
implemented by a so-called two- mode Noon state (Kok et al, 2002): 

|^Ar) = i=(|7V,0) + |e^^'^0,7V)) , (18) 

where \k) is an fc-photon state. When we choose Eat to be 

J:n = \N,0){0,N\ + \0,N){N,0\, (19) 
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then the uncertainty in the phase is again given by A(p — 1/N. However, there is also 
something else going on . . . 

When we calculate the expectation value of Sjv, we find that it varies as cos 1^9. 
That is, the geometric distance between the two points that can be distinguished 
according to the Rayleigh limit is now N times smaller than in the classical case! As 
a consequence, we can in principle read and write much smaller features with this 
technique. This is called quantum lithography (Boto et ai, 2000), and we will now 
give a full description of this phenomenon in terms of the quantum interference of two 
optical modes. 



3.1. Interference on a surface and the Rayleigh limit 

Suppose two plane waves characterised by fci and fc2 hit a surface at an angle 6 from 
the normal vector. The wave vectors are given by 

fci = fc(cos 0, sin 0) and k2 = k{cos9, — sm9) , (20) 

where we assume |fci| = |fc2| = k. The wave number k is related to the wavelength of 
the light according to fc = 2tt/X. 

In order to find the interference pattern in the intensity /, we sum the two plane 
waves at position r at the amplitude level: 



I{r) oc 



When we calculate the inner product (fci 
expression 

I{x) oc cos^(fca;sin6') 



4cos^ 
fca) • r 



(fci - ^2) • r 



(21) 



'/2 from Eq. we obtain the 

(22) 



for the intensity along the substrate in direction x. 

As we saw above, the Rayleigh limit is given by the minimal resolvable feature size 
Aa; that corresponds to the distance between an intensity maximum and an adjacent 
minimum. From Eq. ^ we obtain 



TT 

kAx sin 9 = — . 

2 

This means that the maximum resolution is given by 

TT TT A 



Ax = 



2fcsin6' 2 sin 6*) 4 sin 6* 



(23) 



(24) 



The maximum resolution is therefore proportional to the wavelength and inversely 
proportional to the sine of the angle between the incoming plane waves and the normal. 
The resolution is thus maximal {Ax is minimal) when sinfl — I, ot 9 = 7r/2. This is 
the limit of grazing incidence. The classical diffraction limit is therefore Aa; = A/4. 
Note that this derivation does not use the approximation sm9 ~ 9, which is common 
when considering diffraction phenomena. 



Quantum lithography & parameter estimation 



6 



3.2. Surpassing Rayleigh's diffraction limit 

So how does quantum lithography work? In the Hmit of grazing incidence, we let the 
two counter-propagating light beams a and b be in the combined entangled state of 
N photons 

I^w)a6 = {\N, 0)ab + e^'^'^lO, N)ab) IV2 , (25) 

where if = kx, with k = 2n/X. We define the mode operator e = (a + b)/V2 and its 
adjoint — {d^ + b^)/\/2. The deposition rate A on the substrate is then given by 
(Boto et al, 2002): 

Aat = (V'w|<5Ar|V'Ar) with = "^^j ' (^6) 

i.e., we look at the higher-order moments of the electric field operator on the substrate. 
The deposition rate A scales with the iV*'' power of intensity. Leaving the substrate 
exposed for a time t to the light source will result in an exposure pattern P{(p) = Ajyi- 
After a straightforward calculation wc see that 

Atv oc (1 + cosiV(^) . (27) 

We interpret this as follows. A path-differential phase-shift in light beam b results 
in a displacement x of the interference pattern on the substrate. Using two classical 
waves, a phase-shift of 27r will return the pattern to its original position. However, 
according to Eq. (|27|l . one cycle is completed after a shift of 2'k/N . This means that a 
shift of 27r will displace the pattern N times. In other words, we have N times as many 
maxima in the interference pattern. These will be closely spaced, yielding an effective 
Rayleigh resolution of Ax = X/AN, a factor of N below the classical interferometric 
result of Ax = A/4. 

We have also shown that quantum lithography can be used to create arbitrary 
patterns in one dimension, rather than just closely spaced lines (Kok et al, 2001). 
This can be achieved by using superpositions of the more complicated state 

IV^at™) = -i= (e^'^l - m, m) -t- e^(^-™)'^|m, N ~ m)) . (28) 

There are two fundamentally different ways states of the form in Eq. H28|l can be 
superposed. We can either sum over the photon number A^: 

N 

|*„.) = ^a„|Vnm) , (29) 

Tl = 

with an complex coefficients, or we can sum over the photon distribution m: 

LJV/2J 

I^-JV) = J2 "mlV'JVm) , (30) 

where [A^/2J denotes the largest integer I with I < N/2 and am again the complex 
coefficients. Every branch in this superposition is an TV-photon state. 

These techniques not only allow us to create arbitrary one-dimensional patterns, 
but extensions to four- mode states also facilitate two-dimensional patterns (Kok et al. , 
2001). Furthermore, by choosing the coefficients of Eq. ()30f) in a special way, Bjork 
et al. showed that one can construct a subwavelength resolution pixel state. This state 
can subsequently be used to illuminate a surface in order to etch an arbitrary pattern 
(Bjork et al, 2001). 
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3. 3. Demonstration of quantum lithography and the preparation of states 

The principle of quantum lithography has been demonstrated using two-photon path- 
entangled states generated by parametric down-conversion (d'Angelo et al., 2001). 
In this experiment, the photon pairs are created in a BBO crystal (/3-BaB2 04). 
Immediately behind the crystal, a double-slit aperture is placed, which blocks most 
pairs. However, when two photons do get through, they are extremely likely to have 
passed through the same sht. 

This way, the state of the light field just after the slits is approximately 
|2, 0)ab + |0,2)^B, where A and B denote the two slits. This state is of the form 
of Eq. H25() . and can therefore be used to beat the Rayleigh limit. By scanning the 
output field with a set of two detectors, the diffraction pattern conditioned on a two- 
fold detector coincidence was mapped out. This pattern was twice as narrow as the 
single-photon diffraction pattern, and the principle of sub-Rayleigh resolution due to 
two-photon quantum lithography was thereby demonstrated. Subsequent experiments 
with photon pairs have confirmed these findings (Edamatsu et al., 2002; Shimizu et 
al, 2002). 

However, in order for quantum lithography to work, we also need a suitable photo- 
resist. In other words, we need a material that is sensitive to multi-photon events in 
order to etch images on its surface. Last year, Korobkin and Yablonovich exposed 
commercial photographic film to photon pairs, and produced the coveted two-fold 
resolution enhancement (Korobkin and Yablonovitch, 2002). Since this experiment 
includes both the state preparation and the imaging process, it qualifies as the first 
complete demonstration of quantum lithography. 

Shortly after the experiment of d'Angelo et al., Cataliotti et al, performed 
frequency measurements using multiphoton Raman transitions in Rubidium atoms 
that are confined in an optical dipole trap (Cataliotti et al, 2001). The multiphoton 
events of up to fifty photons resulted in a spectral width that is below the Fourier 
limit. Although not exactly quantum lithography, this experiment strongly suggests 
that the quantum Rayleigh limit of X/AN is correct. 

Another essential ingredient for quantum lithography is the generation of the 
required quantum states of the light field. The states in Eqs. (j28|l and H29|l are very 
complicated, and it is not quite clear how they can be generated efficiently without 
large Kerr nonlinearities (Gerry and Campos, 2001). The production of A^-photon 
entangled states conditioned on non-detection was proposed (Fiurasek, 2002), as well 
as the creation of Noon-states based on single-photon detection (Lee et al, 2002; Kok 
et al, 2002; Gerry et al, 2002). 

4. Conclusions 

We have shown that quantum lithography and Heisenberg-limited parameter 
estimation are two manifestations of the same principle: Instead of many separate 
measurements of ip that lead to the shot-noise limit (in, for example, an experiment 
using N trials with the single-photon path-entangled state |1, 0) +e'^'^\0, 1)), entanghng 
the resources of these N measurements to conduct a single-shot experiment (e.g., using 
the A^-photon Noon-state | A^, 0) -I- 10, A)) can reduce the noise to the Heisenberg limit. 
Similarly, instead of Rayleigh limited single-photon diffraction patterns, we can use 
A^-photon entangled states to increase the resolution by a factor A^. Mathematically, 
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the hyperresolution of quantum lithography appears naturally in the derivation of 
Heisenberg-limited parameter estimation. 
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